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Weyl semimetals are three dimensional analogs of graphene where the energy of the excitations are 
a linear function of their momentum. Pyrochlore Iridates {A-zIriOi) with A =yttrium or lanthanide 
element) are conjectured to be examples of such a system, with the low energy physics described 
by twenty four Weyl nodes. An intriguing possibility is that these materials provide a physical real- 
ization of the Adler-Bell-Jackiw anomaly. In this letter we investigate the properties of pyrochlore 
iridates in an applied magnetic field. We find that the dispersion of the lowest landau level de- 
pends on the direction of the applied magnetic field. As a consequence the magneto-conductivity in 
an electric field, applied parallel to the magnetic field is highly anisotropic, providing a detectable 
signature of the semi-metallic state. 



Graphene [l[ and topological insulators 043| have pro- 
vided the venue for condensed matter realizations of non- 
trivial phenomena originally in the realm of high energy 
physics. Massless relativisitic fermions|5J, Klein tun- 
neling Q, Theta vacuum (i.e. axion electrodynamics) 
and Majorana modes [a, 0] are a few examples. A 
common feature of these systems is that the low energy 
physics is described by a two component Dirac Hamil- 
tonian with the fermionic momentum is confined to two 
dimensions. Recently pyrochlore iridates have been con- 
jectured to realize the Adler-Bell-Jackiw (ABJ)d, Q chi- 
ral anomaly, adding another example to the growing list. 

Wan et al. [l(J explored the possibility of the three 
dimensional analog of graphene being realized in the py- 
rochlore iridates. These materials have a large spin or- 
bit couplings and are in a regime of intermediate cor- 
relations, making them promising candidates to realize 
topological insulators [ll[ [l2[ . They have a magnetic 
ground statefH, [3 and, within a LSDA+U+SO calcu- 
lations, conjectured to be semi-metals. Most strikingly 
the low energy physics is described by the Weyl equa- 
tion, which is the two component version of the Dirac 
equation. There are 24 Weyl nodes, three around each 
L point ([111] and equivalent directions) in the Brilloiun 
zone (see figH]). Nodes related either by inversion or re- 
flection about the {xy, yz, zx} planes have opposite chi- 
rality. Consequently the material is expected to have 
an anomalous hall response to applied uniaxial pressure 
and is susceptible to charge ordering in large magnetic 
fields [Hj]. 

In a quantizing magnetic field, the Lowest Landua 
Level (LLL) is a linear function of the magnitude of a mo- 
mentum, with the sign determined by the band structure. 
If an electric field is applied parallel to the magnetic field, 
the Adler-Bell-Jackiw (ABJ)[|[g| axial anomaly leads to 
an anomalous magneto-conductance. The origin of this 
effect is in the production of Weyl fermions of a given 
chirality and an equivalent annihilation of the opposite 
chirality [liij]. This translates to a transfer of particles 
from one Weyl node to another of opposite chirality at a 
constant rate. To reach a steady state, this is balanced 
by inter-node scattering due to impurities. 




FIG. 1: The low energy physics of pyrochlore iridates is de- 
scribed by linearly dispersing fermionic modes near nodes in 
the band structure. The position of Weyl nodes in the Bril- 
louin zone is shown in the figure ( figure taken form ref[10|]). 
There are three nodes, with the same chirality, located in the 
vicinity of the L points, and nodes related by inversion have 
opposite chirality. 



In the absence of the magnetic field, the intra-nodc 
scattering is quite effective in relaxing the momentum. 
In the presence of a large field, such that only the LLL 
is occupied, intra-node scattering is suppressed due to a 
lack of phase space. The only scattering mechanism avail- 
able are processes that involve different nodes of opposite 
chirality. Since these nodes are located at different point 
in the Brillouin zone, Nielsen and Ninomeya (l6| argued 
that the corresponding scattering rate is much smaller 
implying a large magneto-conductivity. 

In this letter we focus on the anomalous magneto- 
conductivity expected in Weyl semi-metals 
general the Hamiltonian takes the form ±<f-V-cr where q is 
the momentum, V is a real matrix, ± label right handed 
(RH) and left handed (LH) chirality, and a = {cr x , a y , 
a z }. We first address the question of the dispersion of 
the LLL for arbitrary V. We find that the energy is de- 
termined by the component of the momentum parallel to 
the applied field. Crucially the velocity of the mode also 
depends on the direction of the magnetic field leading to 
highly anisotropic response. 

The low energy Hamiltonian in the vicinity the Weyl 
nodes of the proposed topological semimetal state is 
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H(q) = A 



(1) 



2mi 2m 2 
+ (j3q z + ciq\ cos (36)) <J y + c 2 q\ sin (30) a x 

where the local z— axis is taken along the T — L direction 
and the local x— axis along the L — K direction in the 
Brillouin zone. Notice that these direction rotate with 
respect to the global coordinates from one L point to an- 
other. At the Weyl nodes, there is a degeneracy among 
two states with opposite symmetry under inversion. The 
physics of the two state system is captured by the a ma- 
trices. Define qo z and <fo_i_ as the displacements of the 
node away from the L point parallel and perpendicular 
to the r — L line respectively. 9 is defined as the angle 
<fo_L makes with the x— axis. The nodes are located at 
9 = pn/3, qQ Z = ±CiqQj_/(3 (positive (negative) for odd 
(even) values of p ) , and qo± satisfying the equation A 
+q% z /2m 1 -qlJlmi = 0. 

In the vicinity of these nodes the Hamiltonian can be 
expanded and written as 



Hi (5q) ^Sq-Vt-a 



(2) 



where Sq = q — qo, and V,; is a real matrix whose en- 
tries depend on the node index i. For the Hamiltonian in 
eqll]the matrix V is generically not symmetric and does 
not have all real eigenvalues. This is quite unlike the 
case considered in the context of the ABJ anomaly [lij 
or the anomalous hall effect [l5(. For a diagonal V ma- 
trix, the direction of the momentum, whose magnitude 
determines the energy if the lowest Landau level, is par- 
allel to the applied magnetic field. Since momentum does 
not commute with the electromagnetic vector potential, 
and the fact that the latter is always transverse to the 
B field, suggest that this property is generally valid. To 
understand the nature of the ABJ anomaly in iridates, 
we first verify that this conjecture holds for arbitrary V 
matrices. 

The procedure for computing the energy is sketched 
out here (details provided elsewhere [17[). Consider 
a single Weyl node in an applied magnetic field, 
B. We can always rotate the coordinate system so 
that the B lies in the xz— plane. In this reference 
frame B = B {sin (9) , 0, cos (9)}. Using a Landau 
gauge we write the corresponding vector potential as 
B {— cos (9) y, 0, sin (9) y}. Given this choice, the sys- 
tem is translationally invariant in the x and z direc- 
tions. The wave-function of the LLL has the form 
{u,v} 4>(y) e~ lSkxX ~ l5kzZ , where u and v are constants. 
Rather remarkable, the energy eo, calculated with this 
choice of gauge, can be written in a gauge invariant form 
as 



eo 



Det [V] 



\\adj [V] • B\ 



-Sq- B 



(3) 



where Det [V] and adj [V] are the determinant and ad- 
jugate of the matrix [V]. The sign of the dispersion, and 
hence the chirality, is determined by the determinant. 
The energy is inversely proportional to the projection of 
the deviation of the momentum in the direction of the 
applied field. 

An important point to note is that the energy depends 
on the full matrix V and not just its determinant. Ap- 
plying to graphene in a perpendicular magnetic field, it 
is easy to see that the energy is identically zero. The mo- 
mentum is in-plane and has no projection along the ap- 
plied field. For isotropic Weyl fermions where V = ±u/I, 
I being the identity matrix, the dispersion reduces to 
±vp6q ■ B. If V is diagonal but anisotropic, the units 
of length can be rescaled to bring the system back to an 
isotropic form [15j . In general the only statement that 
can be made is that the mode disperses only in the direc- 
tion of the applied field, while the actual velocity cannot 
be determined without a detailed evaluation of the de- 
nominator in eqnJ3J 

Let us now consider the case of the pyrochlore iridates. 
Given the Hamiltonian (eqnfl]) we can construct the rel- 
evant matrices. Since the different Weyl nodes in the 
iridates have different V, all constructed with respect to 
their own local coordinate system, a general solution is 
prohibitive. To make further progress we use parame- 
ters that best fit the LSDA+U+SO calculations fill : 
mi= m 2 = 0.5 eV~\ a = c 2 = 1.0 eV, (3 = 0.5 eV, A 
= 0.18 eV (q is dimensionless) . These parameters give 
q a ± = 0.48 and q Qz = ±0.22. 

To get insight on the general properties let us look 
at the nodes near [111]. The matrices for 9 = 
{0, 2tt/3, -2tt/3} are 



V fl 



V e=0 = 



=±2tt/3 




(4) 



qo± 

±^q 2 01 _ TV3q 
2<z, 



tu 



Since the system possesses three fold rotation symme- 
try about the T — L axis, these matrices are related by 
120 degree rotation about the z— axis. If U is a rota- 
tion matrix for 27r/3 rotation about the local z— axis, 



than UV 0=O = V 0=2 7r/3, UV e=2 7r/3 = V 6 



-27r/3 



and 



UVg = _27r/3 = Ve = o- In other words, knowing one is 
sufficient to generate the others. The determinant of 
these three matrices are all equal and given by Det[V]i 

The adjugate of a matrix V is the transpose of the ma- 
trix of the cofactors of V. The adjugates corresponding 
to those in eqn[5]are 
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adj[V 8=0 ] = 



adj[V e=±2 T T / 3 



2%± 



qo± + ®qoj_qoz 
3ql±qoz 



-3 2 
4%-L 



I (qo± + 6qlj_qoz 

i 3\/3 2 
^ 4 



(5) 



-6<Z 3 , 



0_L 



The adjugate matrices also respect a rotation proper- to the density of states and is given by 



ties about the z— axis. In this case we find adj [V] e=0 U 



adj [V] 



0=2tt/3' 



6>=4tt/3 



u 



t _ 



dj[V] e=2n/3 Ul = adj[V} e= 
idj [V] 



6»=4tt/3 



and 



Thus, unlike the Hamil- 



tonian, the adjugates satisfy a right multiplication prop- 
erty. 

We can construct all other matrices by symmetry. The 
matrices for the Weyl points related by inversion to those 
near [111] are obtained by changing the sign of the third 
column or equivalently the sign of a z . This remains true 
for all nodes related by inversion as the local z— axis 
changes sign. The nodes near [111] have the same struc- 
ture as those near [111] while those near [111] have the 
sign of the first column changed. This ensures the geom- 
etry and helicity obtained within LSDA+U+SO calcula- 
tions [l(| ■ Given that the adjugate matrices are different 
at different nodes, we expect the response to a magnetic 
field to be highly anisotropic. To explore this further we 
turn to the magneto- conductivity. 

Since the dispersion of the LLL is in the direction of the 
applied magnetic field, only the response to an electric 
field, E = EqB applied parallel to the magnetic field 
will be considered. We will assume that the magnetic 
field is strong enough so that only the LLL is occupied. 
The dispersion being linear, elastic scattering within a 
single node is suppressed due to the lack of phase space. 
Stated differently no momentum relaxation is possible 
within a node. However the scattering between two nodes 
cannot be ignored. In fact, scattering between cones of 
opposite chirality is necessitated by the ABJ anomaly. 
In the presence of the electric field, there is a generation 
of RH particles and annihilation of LH particles. The 
rate of production is given by jV| | = 5qeB/h = e 2 EoB/h, 
where eB/h is the degeneracy of the LLL and 5q\\ is the 
component of the momentum parallel to the magnetic 
field. This rate has to be balanced by the scattering 
between nodes with the opposite chirality to maintain a 
steady state. 

Weyl nodes always come in pairs with opposite chiral- 
ity. Before considering all 24 Weyl nodes, we first look 
at a single pair whose energies are given by eo = ±v5qu. 
For impurity scatterers, where the transition probability 
between states with different moment is independent of 
their momenta, the scattering rate, r" 1 is proportional 



l_ = c e_B 
t vn 



(6) 



where C is a constant determined by the strength of scat- 
tering potential. We will comment on more general forms 
of scattering later. The ABJ anomaly requires an imbal- 
ance in particle number at the RH and LH Weyl nodes. 
If the difference in chemical potential between the RH 
and LH nodes is A/z, than energy balance requires 



A/i = cEqvt 



(7) 



This is because over the scattering time r, the momentum 
changes by eE$T. Since the change in energy is v times 
change in momentum, the net energy transferred from 
one cone to the other is cEqtv. A steady state is achieved 
if this transfer is balanced by the difference in chemical 
potential. 

The difference in chemical potential leads to a current 
given by 



J A = e(eB/h)An 
— e(eB/h)eEovr 
= C- 1 e 2 v 2 E 



(8) 



The subscript A refers to the anomalous response. So far 
we have only considered momentum independent scatter- 
ing rate. Let us look at the response in the presence of 
screened charged impurities characterized by matrix ele- 
ments of the form C/(q 2 + k 2 ). The internode scattering 
probes intermediate to large momenta, as the nodes are 
physically separated in the Brillouin zone. The transi- 
tion probabilities, which are proportional to the square 
of the matrix elements, fall of as C/|<7| 4 giving a large 
t. In contrast, the conductivity in zero magnetic field is 
much smaller as it is dominated by intra-node scatter- 
ing (~ C/k 4 ). Thus the conductivity in the presence of 
the magnetic field can be much larger than that in zero 
field. The results in eqlUand eqEl and the argument of 
large magnetoconductance, are the main conclusions of 
Nielsen and Ninomeya [16[ . 
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Having established the expected nature of the 
magneto-conductance, we return to the discussion of the 
iridatcs. The key feature of eqJ5]is its dependence on the 
velocity v. As we have seen in our discussion of the dis- 
persion of the LLL, the velocity can be tuned by changing 
the direction of the applied field. Moreover the different 
Wcyl nodes have different Fermi wave-vector and differ- 
ent density of states. We will exploit the fact that the 
nodes always come in pairs and generalize the formulas 
derived above. We focus on three orientations of the ap- 
plied field and assume that the dominant scattering is 
among nodes related by inversion. 

(1) Magnetic field along a principle axis (say z): All the 
nodes on the Y — L — K plane have the same dispersion of 
the LLL. For the parameters specified above, the energy 
is £o ~ ±0.14<5p • B, where the sign is determined by the 
chirality. All other nodes have eo s» ±0.38<5p • B. The 
velocity of low energy fermions differ by almost a factor 
of 3. Labeling the velocities as v\ = 0.14 and V2 — 0.38 
the anomalous current generated is 

J\ = C~ x e 2 (8vl + 16w|) E (9) 
= 2A7C- 1 e 2 E 

(2) Magnetic field along Y — L (say [111]): For this case 
the nodes divide into three categories. The six Weyl 
points that are in the vicinity of [111] have the same 
velocity given by v\ — ±0.15. At the other six equiva- 
lent L points, the three Weyl nodes split into one with 
velocity V2 = ±0.68 and two with velocity V3 = ±0.24. 
The larger value is associated with the node whose go's 
have a positive projection along the magnetic field. The 
anomalous current is 

J| = C*-l e 2 ( 6w 2 + Q V 2 + 12u 2) Eq (10) 

= 3.60C -1 e 2 .Eo 

(3) Magnetic field along Y — K (say [110] J: The nodes 
again split into three groups. There are four with velocity 
V\ = 0.69 corresponding to the nodes which live in the 
r — L — K plane perpendicular to the magnetic field. 
The other eight nodes close to these four have the same 
velocity of V2 — 0.3. The four nodes in the r — L — K 
plane containing the magnetic field have V2 = 0.3 and 
the remaining eight have i>3 = 0.15. 

j3 = c*-! e 2 (Avl + 12*| + 8W3) £ (11) 
= 3.16C _1 e 2 .E 

It is clear from eqns. l9l ITOl and [TT1 the proposed band 
structure of pyrochlore iridates leads to an anisotropic 
magnto-conductivity. In particular the anisotropy, de- 
fined as the difference in conductivity divided by the sum 
( J\ —J a) I \ J\ +J\) 1 i s as l ar g e as 19% between a field 
along a principal axis and one along Y — L. 



In analyzing the results of the three cases considered, 
we find that the density of states are almost equal. Thus 
thermodynamic properties such as specific heat are not 
as sensitive to the dependence of the dispersion on the 
magnetic field (in particular its direction). The fact that 
the anomalous production of Weyl fermions leads to a 
current that depends on the velocity provides a clear sig- 
nature of the novel semi-metallic state. If one relaxes 
the assumption of scattering only between nodes related 
by inversion, than r is a constant since the density of 
states does not change. In this case the current is pro- 
portional to sum of the absolute velocities of the nodes 
and the anisotropy is reduced to 5%. The presence of 
topologically protected surface states connect bulk nodes 
in pairs providing a mechanism for the transfer of Weyl 
particles [10j. How the surface states and the nature of 
scattering affect the anisotropy will be the subject of a 
future study. 

In this letter we have considered the nature of 
magneto-conducatnee in Weyl semi-metals. In the pres- 
ence of an applied field the energy of the LLL is obtained. 
We find that the the level always disperses linearly with 
respect to the momenta parallel to the applied field. This 
result is general and holds even for systems with arbitrary 
Weyl hamiltonians. On the other hand the velocity of the 
low energy modes depends on the direction of the applied 
field. This feature is exploited in the context of the py- 
rochlore iridates. Using the symmetries of the crystal, we 
derive the Weyl equation at the 24 nodes. Applying mag- 
netic fields in different directions allows us to manipulate 
the low energy physics, leading to anisotropic transport 
properties. 
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